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ABSTRACT 

We present the first estimate of age, stellar metallicity and chemical abundance ratios, 
for an individual early-type galaxy at high-redshift {z = 1.426) in the COSMOS field. 
Our analysis is based on observations obtained with the X-Shooter instrument at 
the VLT, which cover the visual and near infrared spectrum at high (R > 5000) 
spectral resolution. We measure the values of several spectral absorptions tracing 
chemical species, in particular Magnesium and Iron, besides determining the age- 
sensitive D4000 break. We compare the measured indices to stellar population models, 
finding good agreement. We find that our target is an old (t > 3 Gyr), high-metallicity 
([Z/H] > 0.5) galaxy which formed its stars at Zform > 5 within a short time scale 
~ 0.1 Gyr, as testified by the strong [a/Fe] ratio (> 0.4), and has passively evolved in 
the first >3 — 4 Gyr of its life. We have verified that this result is robust against the 
choice and number of fitted spectral features, and stellar population model. The result 
of an old age and high-metallicity has important implications for galaxy formation and 
evolution confirming an early and rapid formation of the most massive galaxies in the 
Universe. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: high-redshift - galaxies: 
stellar content. 


1 INTRODUCTION 

The cosmic history of galaxy mass assembly represents one 
of the open key questions in cosmology. Early-type galax¬ 
ies (ETGs) are the most effective probes to investigate this 
topic, as they are the most massive and oldest galaxies in the 
local Universe and most likely those whose stars formed ear¬ 
liest. Observations have shown that a population of massive 
and passive galaxies is already in place a t high redshift, when 
the Universe was o nly a few Gyr old (iGimatti et aLll2004l . 
ISaracco et ahllioOSl) . So far, the main physical parameters 


* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile 
(program: 086.A-0088(A)) 
t E-mail: ilaria.lonoce@brera.inaf.it (OAB) 


related to their formation and assembly have been mainly 
estimated on local ETGs, and their ageing and evolution can 
mix up and confuse the original properties when the bulk of 
their mass formed and assembled. Information on the star 
formation (SF) time-scale of high-t; ETGs can be obtained 
from the det ailed chemical abund ance ratios of their stellar 
populations (iThomas et al.ll2005l) . which can be derived by 
a detailed spectral analysis. Indeed, the abundance of Iron 
with respect to a-elements is tightly correlated with the time 
delay between Type I and Type II supernovae (SN), giving a 
direct probe of the time-scale within which SF has occurred. 

Up to now, only few works have e xperimented a spec¬ 
tral analysis on ETGs at z > 1 (iLonoce et al.l l2014l . 

lOnodera et al.l l2012l . Ijorgensen et al.l l2014l) due to the low 
S/N of the available spectroscopic data, and they were 
mostly focused on age estimates, in particular using the UV 
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region dCimatti et alJlioOSh . The analysis of the rest-frame 
optical spectrnm is still lacking. 

Furthermor e, measnres are n sually performed on 
stacked spectra llOnodera et al.l[2015l l. thus deleting possible 
peculiarities of single objects. 

A single-object measurement of age, stellar metallicity 
and chemical abundance ratios of 2 > 1.2 ETGs is miss¬ 
ing at the present time. We fill this gap, presenting the hrst 
attempt to measure the detailed chemical composition, be¬ 
sides age, of a 2 ~ 1.4 ETG directly in the early stages of 
its evolution. 

Throughout this paper, we assume a standard cosmol¬ 
ogy with Ho = 70 km s“^ Mpc“^, Qm = 0.3 and Qa = 0.7. 


2 COSMOS-307881: SPECTROSCOPIC DATA 


Our target is a bright and massive (> 10^^ M©) ETG 
from the K-selected galaxy catalo g in the GOSMOS (Cos - 
mological Evolution Survey) field ( McCracken et aLll2010ll . 
It is one of the 12 galaxies with Ks(Vega) < 17.7 se¬ 
lected by iMancini et al.l (l2010ll on the basis of three cri- 
teria: (^i) non-detecti on at 24 ^m in the S'pjt 2 er-|-MIPS data 
dSanders et al.l [20071 1: (ii) visual elliptical morphology (see 
Fig. [T1 top panel); (iii) multicolor Spectral Energy Distri¬ 
bution (SED) consistent with old and passive stellar pop¬ 
ulations with no dust reddening. All availab le information 
are sh own in Tabled] For furthe r details see IMancini et aP 
(l2010li and lOnodera et al] (l2012tl . From these studies we can 
infer that our target is an old and slightly compact (Re ~ 0.3 
Rr°, w.r.t. the local size-mass relation) ETG, and most of 
its mass is composed by a passively evolving stellar popula¬ 
tion. In the following we present a new analysis of the stellar 
population properties of this ETG based on new spectro¬ 
scopic data, which, thanks to the high resolution in a wide 
spectral range of the X-Shooter spectrograph, allows us to 
measure the stellar metallicity of this high -2 target. 

Spectroscopic observations were car ried out with the X - 
Shooter spectrograph on the VLT/UT2 llVernet et al.l[201ll l 
in Italian guaranteed time during the nights 9 — 10 Febru¬ 
ary 2011, (program: 086.A-0088(A)). The target has been 
observed under bad sky condition (seeing > 2 ") during the 
first night (about 1.7h for the VIS and UVB arms and about 
1.9h for the NIR arm), while during the second night more 
than 4 hours of exposure have been collected under good 
sky conditions (seeing ~ 0.8”). We decided to consider only 
the latter set of observations. The use of the l”.0 slit in the 
UVB arm and 0”.9 slit in the VIS and NIR arms, resulted 
in a spectral resolution of 5100, 8800 and 5600 in the UV, 
VIS and NIR respectively. 

The data reduction has been performed taking advan¬ 
tage of the ESO pipeline l|GoIdoni et al.lf2006l l regarding the 
first steps of the process, and completed by means of IRAF 
tools. In Fig. [T] the reduced mono-dimensional spectrum of 
307881 (black line) is shown in compa rison with the pho¬ 
tom etric points (cyan diam onds, from l^uzzin et ^l2013l . 
and lMcCracken et al.l[201^ '). Overlaid is a 4 Gyr, super-solar 
metallicity (Z— 0.04) simple stellar population model of 
iMaraston fc Str ombackI (|2(ill|) (Mil) based on the MILES 
stellar library ( Sanchez-Blazauez et (red line). Note 

that this template is the most similar in term of age and to¬ 
tal metallicity, to the best-fit solution that will be derived 


in the next section. The full spectr al model, as known, does 
not in clude the a/Fe parameter ('see lMaraston fc StrombackI 
(l 201 lll '). which we shall derive using selected a b sorpt ion-line 
models in the next section fsee lThomas et aP (1201 lI U. 

The original high resolution of the spectrum (e.g. ~ 1 A 
in the VIS) has been decreased to match the lower resolution 
of the models (i.e. ~ 2.54 A at rest-frame, and ~ 6 A in the 
VIS at 2 = 1.4), which are used to perform the analysis of 
spectral indices (for details see Lonoce et al. 2015 in prepara¬ 
tion). Note that we adopted the MILES-based Mil models 
instead of the higher resolution ones (Mll-ELODIE or Mll- 
MARCS, see www.maraston.eu/Mll) to take advantage of 
the increasing S/N in downgrading the observed resolution. 
The value of the signal to noise ratio (S/N) obtained in the 
spectral region around 5000 A restframe is ~ 7 per pixel. 

In Fig. Owe show how the H/3 and Mg 6 are expected 
to appear on a synthetic template (upper first panel) that 
broadly reproduces the observed one. The model is firstly 
adapted to the measured velocity dispersion (upper second 
panel), then downgraded to be as noisy as the observed spec¬ 
trum (upper third panel) and finally compared to the ob¬ 
served one (bottom panel). The similarity between expected 
and observed absorption features is quite evident. 


3 SPECTRAL ANALYSIS 

As a first step, we measured the redshift of 307881 fitting 
the Mgt line region, which is the cleanest from the back¬ 
ground residuals, as it can be seen in Fig. [T] (bottom panel), 
finding 2 = 1.426 ±0.001. The Mgt, line region has been also 
used to find a best fitting velocity dispersion estimate, that 
resulted to be a = 385 ±85 km/s. More details on the fitting 
procedure adopted to fix both 2 and a will be described in 
the forthcoming paper. 

We then selected some indices whose absorption fea¬ 
tures are clearly visible in the observed spectrum, to try to 
simultaneously derive both the mean age and the metal¬ 
licity of its s tellar populatio n. The selected indices are: 
D4000 index (iHamiltonl 1 1985l l . H 7 F (IWorthev fc Ottavianil 
Il997fl . G4300, Fe43 83, Ga4455, Fe4531, H/3, Fe5015 and Mg 6 
(Lick/IDS system. IWorthev et al.l[l994l L In particular, it is 
well known that the Mgt index is the best metallicity and 
chemical abundanc e dependent inde x in the region around 
5000 A restframe jKorn et al.l 1200^ 1. In Table [2] we report 
the measured values of these indices together with their er¬ 
rors derived by means of Monte Carlo simulations set on the 
uncertainties in the flux measurements. 

Finally, we want to point out the presence of the 
[OH] 3727 emission line (Fig. [T] middle panel). We can ex¬ 
clude that its origin is due to an active AGN, since we do 
not see any other signature in the observed wide spectral 
window. There are reasonable possibilities that this emis¬ 
sion is caused by the UV i onizing emission of old stars in 
post main-sequence phases jYi fc Yoonll2004h . as confirmed 
by UV indices (Lonoce et al. 2015 in preparation), while a 
strong contribution from SF can be excluded. Indeed, as it 
can be noted from Fig. 0 (left panels), it is highly unlikely 
that the H/3 feature is affected by emission, considering also 
that its value suggests a stellar population age in good agree¬ 
ment with that derived by the H 7 and D4000 index. 
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Figure 1. COSMOS 307881. The galaxy observed spectrum (black line) is compared to a model llMaraston &: Strombacj|201lh with 
age of 4 Gyr and super-solar metallicity (Z= 0.04), and to observed photometric data (cyan diamond). Top panel: VIS and NIR spectral 
region together with the HST/ACS I-band image of the target. Middle panel: zoom of the 4000 A rest-frame region. Bottom panel: zoom 
of the 5000 A restframe region. The main absorption (and one emission) lines in each spectral region are highlighted. Dark green lines 
indicate the residual spectrum. 
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Table 1. COSMOS-307881. Data derived from the analysis of lOnodera et ahl l|2Q12l~) : K-band magnitude in Vega system (Ks); 
spectroscopic redshift s tellar popu lation age (Agep^ot) and logarithm of the stellar mass (logAf*) derived from SED 

fitting assuming a Chabrier IMF llChabrierll2003h : effective radius (Re); degree of compactness (C=Re/Re,z=o); Sersic index 
(n). Units of right ascension are hour, minutes and seconds, and unites of declination are degrees, arcminutes and arcseconds. 


ID 

RA 

DEC 

K, 

H^phot 

logA4* 

Re 

C Sersic n 




(Vega) 

(Gyr) 

(Mo) 

(kpc) 


307881 

10:02:35.64 

02:09:14.36 

17.59 

1.4290±0.0009 3.50 

11.50 

2.68±0.12 

0.32 2.29±0.10 


(*)This work: Zspec = 1.426 ± 0.001 



1170 1175 1180 1185 

Wavelength [nm] 



Figure 2. H/3 (left panels) and Mg;, (right panels) features in the 4 Gyr, 2ZO metallicity model shown in Figure 1. From top to bottom 
we show: model at 2 ~ 1.4; model corrected for a = 385 km/s; model downgraded for the observed poissonian noise; model compared to 
the observed spectrum (point-dashed red line). 


Table 2. Measured indices values. 


4 MODEL COMPARISON 


Index 

value 

D4000 

2.44 ± 0.12 

D„4000 

2.42 ± 0.17 

H7f 

-1.56 ±0.92 

G4300 

6.52 ± 1.12 

Fe4383 

7.40 ± 1.74 

Ca4455 

1.06 ± 0.83 

Fe4531 

3.20 ± 1.40 

H/3 

2.52 ± 0.93 

Fe5015 

3.89 ± 1.91 

Mg6 

5.75 ±0.81 


In Fig. |3] we show the observed Mgt and H/3 indices 
(b lue diamond), compa red to the predictions of the SSP 
of iThomas et al.l (I2OIII ) , based on the MILES library for 
a wide range of ages (from 0.1 Gyr shown up to the age of 
the Universe at z ~ 1.4, i.e. 4 Gyr); super-solar metallicities 
[Z/H]= 0.35 (red lines) and [Z/H]= 0.67 (cyan lines), and 
various a/Fe paramete r s, nam ely [a/Fe]= 0.0, 0.3, 0.5. Mod¬ 
els assume a ISalpet^ (Il955l ) initial mass function (IMF), 
and are corrected for the measured velocity dispersion value 
of a — 385 ± 85 km/s. As it can be seen, the extreme value 
of the Mgi, index fully requires high metallicity models up to 
[Z/H]= 0.67. In particular, in Fig. [3]the model expectations 
for these indices are reported also in case of non-solar values 
of the Q-enhancement [a/Fe], from 0.3 to 0.5 (dashed and 
dotted lines respectively). 
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Figure 3. H/3 versus Mg;, plot. Comparis on between th e mea - 
sured indices (blue diamond) and models of iThomas et aO ll201lh 
(lines). Ages run from 0.1 to 4 Gyr, for super solar metallicities 
of [Z/H]= 0.35, 0.67 (red, cyan lines), and a-element abundances: 
[a/Fe]= 0.0, 0.3, 0.5 (solid, dashed, dotted lines). The models are 
corrected for the measured value of ct = 385 km/s. 


Models corresponding to extreme values of [a/Fe]> 0.5 
seem to be required to match the observations. 

The behaviour of the Mg(, index, requiring such extreme 
values of Z, is confirmed also when considering the other 
absorption lines that we were able to measure on this X- 
Shooter spectrum. Indeed, as it can be noticed in Fig. 0] 
where we propose 3 examples of Lick indices (G3400, H/3 
and Fe4383, upper, middle and bottom panel respectively) 
as a function of D4000, all measured indices consistently 
point towards a very-high metallicity ([Z/H] ~ 0.67, cyan 
lines) being only marginally consistent with the 2Zq values 
([Z/H] ~ 0.35, red lines). 

Notice t h at th e highest metallicity models in 
iThomas et al.l (l201ll ) are partly in extrapolation, as 
they are sampling the edge of the parameter space in terms 
of the available empirica l fitting functions for su ch extreme 
stellar parameters (see Ijohansson et al.l il2010l )'). At the 
same time, the underlying stella r trac ks are based on real 
calculations (see iMaraston et al.l (l2003l ) for details). 

More quantitatively, we have computed the best fit¬ 
ting solution obtained comparing all 9 observed indices val¬ 
ues (Table [2]) with models. The free parameters were age 
(0.1 — 4.5, truncated at the age of Universe, with step 
0.1 GyrfI, the total metallicity (from [Z/H]= —2.25 to 
[Z/H]= 0.67, with step 0.01) and the a/Fe-enhancement 
(from [Q/Fe]= —0.3 to [a/Fe]= 0.5, with step 0.01). The 
minimum value corresponds to an age of 4.0/tQ g Gyr, 
metallicity [Z/H]= 0.611 q;q| and [Q:/Fe]= 0 . 454 q; 59 , with 
= 0.7 and an associated probability 70%. Errors indi¬ 
cate the range values of these parameters over all the solu¬ 
tions associated to probabilities larger than 65%. The dis¬ 
tributions of the 3 fitting parameters, displayed in different 

^ Note that we have performed the analysis using all available 
model ages (up to 15 Gyr) thereby ignoring the age of the Uni¬ 
verse as a constrain. However, as no particular improvement was 
noticed in the derived quantities, we decided to focus on ages 
within the age of the Universe at .2 1.4 



D4000 

Figure 4. Comparison between four measured i ndices (G4300, 
H/3 an d Fe4383 versus D4000) and the models of iThomas et ail 
(l201lh . Ages run from 0.1 — 4 Gyr, metallicities from [Z/H]= 0, 
0.35, 0.67 (black, red, cyan lines), and a-element abundances from 
[a/Fe]= 0, 0.3, 0.5 (solid, dashed, dotted lines). The measured 
values are shown with a blue diamond. Indices values of models 
are corrected for the measured value of a = 385 km/s. 


X^ ranges, are shown in Fig. [5] top panel. A global picture 
of the x^ values can be seen in Figure 0 ] where the minimum 
X^ trends of the 3 parameters of all solutions are shown. It is 
easy to notice that ages < 2 Gyr can be completely excluded 
due to the rapid increasing of their x^ values toward younger 
ages. Instead for ages > 4.5 (limit of the Universe age, not 
shown here) the x^ values remains practically constant. 

Furthermore, we found that models with Z^Zq provide 
a fit of the free parameters with a probability less than 0 . 1 %. 

We also verified the strength of this result by repeating 
the same fitting process selecting smaller and different set 
of indices, finding very similar solutions with respect to the 
previous ones based on the whole set of indices. Two exam¬ 
ples are shown in Fig. [S] (middle and bottom panels): the 
distributions of the fitting solutions are obtained from two 
smaller sets of indices (i) D4000, G4300, Hy, Fe4383, H/3, 
Fe5015 and Mgt and ii) D4000, Hy, H/3 and Mgj) which 
lead to a best-fit solution of i) age= 4.0 Gyr, [Z/H]= 0.61 
and [a/Fe]= 0.44 with x^ = 0'9> ^nd ii) age= 4.0 Gyr, 
[Z/H]= 0.60 and [Q/Fe]= 0.5 with x^ = 0.5, both totally 
consistent with the all-indices one. 

We also evaluated the feasibility of this analysis on 
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Figure 7. Distributions of [Z/H] and [o/Fe] derived for a set 
of 500 mock spectra built on a model spectrum with parameters 
(Table [Sjl as those of the best-fit to the observed spectrum and 
adding the poissonian noise. The two distributions are well peaked 
around the true original values (red vertical lines). 


this low S/N spectrum, in particular for metallicity and a- 
enhancement estimates, by repeating it on a set of 500 mock 
spectra built on a model spectrum with parameters as the 
best-fit to the observed one (cfr. Table (S]) and downgraded 
with the poissonian observed noise. The obtained distribu¬ 
tions of the measured [Z/H] and [a/Fe] are shown in Fig. 
[T] They all peak around the true original values (red verti¬ 
cal lines) demonstrating that within the declared errors the 
obtained values are solid. 

Furthermore we tested if the large error on the velocity 
dispersion could affect the results by performing the same 
analysis assuming cr = 300 km/s. We found the same best-fit 
solution (4.1 Gyr, [Z/H]= 0.6, [a/Fe]= 0.41 with — 0.6). 

Finally, in order to test the model dependence of 
this result, we repeated the same analysis adopting the 
iBruzual fc Chariot! (l2003l ) models (BC03). These models do 
not include the [a/Fe] parameter, hence we shall use them to 
constrain age and total metallicity solely. The BC03 models 
cover a slightly lower Z range with respect to our fiducial 
TMJ models, and are based on different stellar evolution¬ 
ary tracks. The BC03 best-fit corresponds to an age of 4.5 
Gyr and a metallicity of [Z/H]= 0.4, which is the maximum 
available metallicity in these models. Hence our result of a 
high-age and high-Z is not model dependant. 

It is important to note the high metallicity is mainly 
derived as a consequence of the maximum allowed age of 
4.5 Gyr. Should we allow the age to be older than the age 
of the Universe, the metallicity will decrease, as a result of 
age/metallicity degeneracy. Actually, the large error bar of 
both the indices and of a prevents a real precise measure 
of the metallicity, but the peculiarity of such a strong Mgt 
absorption band combined with the narrow range of possible 
ages (due to its redshift), make necessary the assumption of 
a very high value of the stellar metallicity. 


5 DISCUSSION AND CONCLUSIONS 

We have performed a detailed spectroscopic analysis of a 
z=1.426, massive {M* ~ early-type galaxy. We 


Table 3. Results: spectroscopic redshift {zspec) and velocity dis¬ 
persion (ct) obtained from this spectroscopic analysis; best-fit val¬ 
ues of age (Gyr), [Z/H] and [o/Fe]. 


Zspec cr (km/s) Age (Gyr) 

[Z/H] 

[a/Fe] 

1.426 ±0.001 385 ±85 4.o/:o;g 

n fii 

U.D1_o.o5 

U.40_o.i9 


gain strong evidence of a high stellar metallicity and a- 
enhancement, and old (relative to its redshift) age (Table|3]). 
These quantities constrain the past SF history experienced 
by the galaxy. In particular, the [a/Fe] ratio, quantifying the 
time delay between Type II SN events, responsible of the 
production of a-elements, and the Type la SNs, related to 
the formation of Fe-peak elements, allows to determine the 
SF time-scale dThornas et ahlBoOSl) . The high value [a/Fe] 
~ 0.4 obtained for 307881 is a direct signature that its SF 
time-scale must have been shor t. In particular , adop ting the 
simple theoretical modelling of iThomas et al.l ll2005l) . where 
the SF is modelled with a Gaussian function, we obtain a 
time-scale At 0.1 Gyr covering the interval within which 
95% of the stars were formed. Considering also the old age 
of its stellar content, this suggests that 307881 formed the 
bulk of its stars at Zform > 5 within a short time-scale of 
At ~ 0.1 Gyr and then passively evolved over the following 
4 Gyr. 

With the high [a/Fe] value suggesting a short SF time 
scale for 307881, the clear indication for an extremely high 
total metallicity opens new issues on the gas enrichment his¬ 
tory of the Universe. It is worth emphasising that the global 
integrated metallicity of the local ETGs never reaches values 
higher than the 1 — 2 Zp). Indeed , we co nsidered the local 
sample analysed in [Thomas et al.l (l2010l ) and quickly veri¬ 
fied that the metallicity and [a/Fe] values of 307881 together 
with its velocity dispersion estimate, are not included in the 
local distribution of values, even if such high values of Z and 
[a/Fe] are expected for dense ETGs as suggested by these 
scaling relations. Thus, this suggests that this galaxy must 
experience mass accretion events (minor merging) from z 
= 1.4 to 2 : = 0 which will move it on the observed local 
scaling relations towards lower values of both Z and a, di¬ 
luting the extreme metallicity stars and confining them in 
the central part of the galaxy. Indeed, such extreme metal¬ 
licity values are in some cases found in t he inner core of lo¬ 
cal massive ETG s llTh omas et ^ l2005l . iTraaer et al.l l2000l . 
iMartln-Navarro et_alj [2(^ ) which are k nown to show metal¬ 
licity gradients I La Barbera et al.ll2012l) . We calculated that 
adding e.g. 15% of a sub-solar metallicity component, as in 
the case of dwarf galaxies, to 307881 would decrease the mea¬ 
sure of its average metallicity to the local observed values. 
On the other hand, the gas metallicities up to now measured 
in z > 3 star forming galaxies result to be solar or sub-solar, 
and do not match at all the high stellar value measured in 
our target gala xy (iMannucci et al l2009l ) and in the centers 
of local ETGs liSnolaor et al. 2008). It is out of the aims of 
the present paper to suggest a possible explanation of this 
missing detection of high- 2 : high metallicity gas. A possi¬ 
bility is that the lower gas metallicity comes from dilution 
through infalling primeval gas, or selective mass loss of met¬ 
als in galactic winds. At the same time, the estimate of the 
metallicity of the Broad Line Regions in quasars at 2 > 4 
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X" distributions 



Age (Gyr) [Z/H] [o/Fe] 

X" distributions 
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Figure 5. Distributions of the fitting-solutions: age, metallicity and ^-enhancement. Colours and symbols indicate solutions in different 
ranges, e.g. green squares include the most probable solutions. Top panel: distributions obtained from the all-indices analysis; middle 
panel: distributions obtained from a smaller set of indices: D4000, G4300, H 7 , Fe4383, Fe5015 and Mg^; bottom panel: distributions 

obtained from a smaller set of indices: D4000, H 7 , H/3 and Mg^. 



Age (Gyr) [Z/H] [o/Fe] 

Figure 6 . Trends of minimum for ^rid [a/Fe] of all obtained solutions. In particular, ages < 2 Gyr can be completely 

excluded due to the rapid degrade of x^ towards younger ages. 
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reveals gaseou s metallicities even higher than the stellar one 
reported here (ljuarez et al ]|200^ . suggesting the possibility 
that their enriched gas is involved in the initial SF events of 
(at least some) high -2 massive proto-elliptical galaxies. 
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